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Abstract. We present a method for numerical computation of conformal 
mappings from simply or doubly connected domains onto so-called canoni- 
cal domains, which in our case are rectangles or annuli. The method is based 
on conjugate harmonic functions and properties of quadrilaterals. Several nu- 
merical examples are given. 

Keywords, numerical conformal mappings, conformal modulus, quadrilater- 
als, canonical domains. 

2010 MSG Primary 30C30; Secondary 65E05, 31A15, 30C85. 

1. Introduction 

In addition to their theoretical significance in complex analysis, conformal map- 
pings are important in classical engineering applications, such as electrostatics 
and aerodynamics [23], but also in novel areas such as computer graphics and 
computational modeling [U [H] . In this paper we study numerical computation 
of conformal mappings / of a domain f2 C C into C. We assume that the domain 
is bounded and that there are either one or two simple (and non-intersecting) 
boundary curves, i.e., the domain Q is either simply or doubly connected. It 
is usually convenient to map the domains conformally onto canonical domains, 
which are in our case rectangles Rh = {z & C : < Re z < 1, < Im z < h} 
or annuli Ar = {z & C : < \z\ < 1}. While the existence of such conformal 
mappings is expected because of Riemann's mapping theorem, it is usually not 
possible to obtain a formula or other representation for the mapping analytically. 

Several different algorithms for numerical computation of conformal mappings 
have been described in the literature. One popular method involves the Schwarz- 
Christoffel formula, which can also be generalized for doubly connected domains. 
A widely used MATLAB implementation of this method is due to DriscoU [S] 
and a FORTRAN version for the doubly connected case is due to Hu [13]. For 
theoretical background concerning these methods see [H [TUl |23]. In addition. 
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there are several approaches which do not involve the Schwarz-Christoffel for- 
mula, e.g., the Zipper algorithm of Marshall [19l |2Q]. A method involving the 
harmonic conjugate function is presented in [121 PP- 371-374], but this method is 
different from ours as it does not use quadrilaterals. For an overview of numerical 
conformal mappings and moduli of quadrilaterals, see [21]. Historical remarks 
and an outline of development of numerical methods in conformal mappings is 
given in [HinilSg. 

In this paper, we present a new method for constructing numerical conformal 
mappings. The method is based on the harmonic conjugate function and prop- 
erties of quadrilaterals, which together form the foundation of our numerical 
algorithm. The algorithm is based on solving numerically the Laplace equa- 
tion subject to Dirichlet-Neumann mixed-type boundary conditions, which is 
described in [11]. To the best of our knowledge, this is the first attempt to 
construct conformal mappings by using hp-FEM. It should be noted, that the 
presented method is not restricted to polygonal domains, and can be used with 
domains with curvilinear boundary as well. 

The outline of the paper is as follows. First the preliminary concepts are intro- 
duced and then the new algorithm is described in detail. Before the numerical 
examples, the computational complexity and some details of our implementation 
are discussed. The numerical examples are divided into three sections: validation 
against the Schwarz-Christoffel toolbox, the analytic example, simply connected 
domains, and finally ring domains. 

2. Foundations of the Conjugate Function Method 

In this section we introduce the required concepts from function theory, and 
present a proof of a fundamental result leading to a numerical algorithm. 

Definition 2.1. (Modulus of a Quadrilateral) 

A Jordan domain Q in C with marked (positively ordered) points 2:1, ^3, -24 ^ 
dfl is called a quadrilateral, and is denoted hj Q = {Q; zi, Z2, z^, Z4). Then 
there is a canonical conformal map of the quadrilateral Q onto a rectangle Rh = 
{Q'; 1 + ih, ih, 0, 1), with the vertices corresponding, where the quantity h defines 
the modulus of a quadrilateral Q. We write 

M(Q) = h. 

Note that the modulus h is unique. 
Definition 2.2. (Reciprocal Identity) 

It is clear by the geometry [THl p. 15] or [211 PP- 53-54] that the following 
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M(g)M(Q) = 1 



where Q = {Q; Z2, z^, Z4, zi) is called the conjugate quadrilateral of Q. 

For basic properties of modulus of quadrilaterals, we refer the reader to [18] and 
[211 Chapter 2]. 

Remark. The identity ([T]) leads to a method for estimating the numerical accu- 
racy of the modulus. For discussion and several numerical examples, see 

2.1. Dirichlet-Neumann Problem. It is well known that one can express the 
modulus of a quadrilateral Q in terms of the solution of the Dirichlet-Neumann 
mixed boundary value problem [T2l p. 431]. 

Let f2 be a domain in the complex plane whose boundary dQ consists of a finite 
number of regular Jordan curves, so that at every point, except possibly at finitely 
many points of the boundary, a normal is defined. Let dQ = AU B where A, B 
both are unions of regular Jordan arcs such that A fl S is finite. Let ipA-, V's be 
real-valued continuous functions defined on A, B, respectively. Find a function 
u satisfying the following conditions: 

1. u is continuous and differentiable in Q. 

2. u{t) = ijjAit), for all t G A. 

3. If d/dn denotes differentiation in the direction of the exterior normal, then 



The problem associated with the conjugate quadrilateral Q is called the conjugate 
Dirichlet-Neumann problem. 

Let 7j, j = 1, 2, 3, 4 be the arcs of dfl between {zi, Z2) , {z2, Z3) , {z^, z^) , {z^, zi), 
respectively. Suppose that u is the (unique) harmonic solution of the Dirichlet- 
Neumann problem with mixed boundary values of u equal to on 72, equal to 
1 on 74, and du/dn = on 71,73. Then by [H Theorem 4.5] or [2n Theorem 



Suppose that Q is a quadrilateral, and u is the harmonic solution of the Dirichlet- 
Neumann problem and let f be a conjugate harmonic function of u such that 
v(Re Z3, Im 2:3) = 0. Then f = u + iv is an analytic function, and it maps Q onto 
a rectangle Rh such that the image of the points Zi, Z2, Z3, z^ are 1 + ih, ih, 0, 1, 
respectively. Furthermore by Caratheodory's theorem [TH| Theorem 5.1.1], / has 



d_ 
dn 



u{t)=lljB{t), 



for all t G -B. 



2.3.3]: 



(2) 
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a continuous boundary extension which maps the boundary curves 71,72,73,74 
onto the hne segments 71,72,73,74, see Figure [ij 




Figure 1 . Dirichlet-Neumann boundary value problem. Dirichlet 
and Neumann boundary conditions are marked with thin and thick 
lines, respectively. 



Lemma 2.3. Let Q be a quadrilateral with modulus h, and let u he the harmonic 
solution of the Dirichlet-Neumann problem. Suppose that v is the harmonic 
conjugate function of u, with v(Re Z3,lm Z3) =0. If u is the harmonic solution 
of the Dirichlet-Neumann problem associated with the conjugate quadrilateral Q, 
then V = hu. 

Proof. It is clear that u are harmonic. Thus v = hu is harmonic, and v 
and V are both constant on 71,73. By Cauchy-Riemann equations, we obtain 
(Vm, Vv) = 0. We may assume that the gradient of u does not vanish on 72, 74. 
In particular, on 74, we have n = Vu/|Vu|, where n is the exterior normal of 
the boundary. On the other hand, on 72, we have n = — Vu/|Vu|. Therefore, 
we have 

dv 1 

= (Vt;, n) = ±— {Vv, Vu) = 0. 
on |Vu| 

By the definition of m, we get 

dv ^du ^ 
On dn ' 

on 72, 74- Thus v and v satisfy the same boundary conditions on 72, 74- Then by 
([1]) and the uniqueness theorem for harmonic functions [21 p. 166], we conclude 
that V = V. □ 



Suppose that f = u-\-iv, where u and v are as in Lemma 2.3 Then it is easy to 
see that the image of equipotential curves of the functions u and v are parallel 
to the imaginary and the real axis, respectively. 
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Finally, we note that the function / constructed this way is univalent. To see 
this, suppose that / is not univalent. Then there exists points 2:1,-22 G Q and 
zi 7^ Z2 such that f{zi) = f{z2). Thus Re/(2;i) = Re/(z2), so zi and Z2 are on 
the same equipotential curve C of u. Similarly for imaginary part, zi and Z2 are 
on the same equipotential curve C of v. Then by the above fact of equipotential 
curves, it follows that zi = Z2, which is a contradiction. 

2.2. Ring Domains. Let E and F be two disjoint and connected compact sets 
in the extended complex plane Coo = C U {00}. Then one of the sets E,F is 
bounded and without loss of generality we may assume that it is E. Then a set 
R = Coo\{E U F) is connected and is called a ring domain. The capacity of R is 
defined by 



where the infimum is taken over all non-negative, piecewise differentiable func- 
tions u with compact support in RU E such that u = 1 on E. Suppose that a 
function u is defined on R with 1 on and on F. Then if u is harmonic, it 
is unique and it minimizes the above integral. The conformal modulus of a ring 
domain R is defined by M{R) = 27r/capi?. The ring domain R can be mapped 
conformally onto the annulus Ar, where r = M{R). In [3] numerical computation 
of modulus of several ring domains is studied. 

3. Conjugate Function Method 

Our aim is to construct a conformal mapping from a quadrilateral Q = {Q; Zi, Z2, Z3, z^) 
onto a rectangle Rh, where h is the modulus of the quadrilateral Q. Here the 
points Zj will be mapped onto the corners of the rectangle Rh. In the standard 
algorithm the required steps are the following: 

Algorithm 3.1. (Conformal Mapping) 

1. Find a harmonic solution for a Dirichlet- Neumann problem associated with 
a quadrilateral. 

2. Solve the Cauchy-Riemann differential equations in order to obtain an an- 
alytic function that maps our domain of interest onto a rectangle. 

The Dirichlet-Neumann problem can be solved by using any suitable numerical 
method. One could also solve the Cauchy-Riemann equations numerically (see 
e.g. [5]) but it is not necessary. Instead we solve v directly from the conjugate 
problem, which is usually computationally much more efficient, because the mesh 
and the discretized system used in solving the potential function u can be used 
for solving v as well. 
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This new algorithm is as follows: 

Algorithm 3.2. (Conjugate Function Method) 

1. Solve the Dirichlet-Neumann problem to obtain Ui and compute the modulus 
h. 

2. Solve the Dirichlet-Neumann problem associated with Q to obtain U2- 

3. Then f = Ui + ihu2 is the conformal mapping from Q onto Rh such that 
the vertices {zi, Z2, z^, Z4) are mapped onto the corners (1 + ih,ih, 0,1) . 

In the case of ring domains, the construction of the conformal mapping is slightly 
different. The necessary steps are described below and in Figure [2] 

Algorithm 3.3. (Conjugate Function Method for Ring Domains) 

1. Solve the Dirichlet problem to obtain the potential function u and the mod- 
ulus M{R) . 

2. Cut the ring domain through the steepest descent curve which is given by 
the gradient of the potential function u and obtain a quadrilateral where 
the Neumann condition is on the steepest descent curve and the Dirichlet 
boundaries remain as before. 



3. Use the Algorithm 3.2 



Note that the choice of the steepest descent curve is not unique due to the implicit 
orthogonality condition. 

4. Implementation Aspects 

The /ip-FEM implementation we are using is described in detail in [TT]. For 
elliptic problems, the superior accuracy of the higher order methods with rela- 
tively small number of unknowns has to be balanced against the much higher 
integration cost and the cost of evaluating the solution at any given point in the 
domain. It should be emphasized though, that the conjugate function method is 
not dependent on any particular numerical PDE solution technique. Indeed, we 
fully expect that similar results could be obtained with, for instance, fine-tuned 
integral equation solvers. 

In the context of solution of the conjugate pair problems, it is obvious that 
we have to integrate only once. Moreover, the factorization of the resulting 
discretized systems can be, for the most part, used in both problems without 
any extra work. Therefore, although in principle two problems are solved, in 
practice the work is almost proportional to that of one. 

However, the computation of the contour lines necessarily involves a large number 
of evaluations of the solution, that also become more expensive as the order of 
the method increases. 
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(a) Ring domain with 
Dirichlet data 0, and 1, 
on the outer and inner 
boundaries, respectively. 



(b) Ring domain: Solu- 
tion of the Dirichlet prob- 
lem with contour lines. 




(c) Conjugate problem 
for the cut domain with 
new Dirichlet data along 
the both sides of the cut. 




(d) Conjugate problem: 
Solution of the conju- 
gate problem with con- 
tour lines. 




(e) Mapped annulus. 
Figure 2. Conjugate Function Method for Ring Domains. 
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4.1. hp -FEM. Here we give a short overview of the /ip-FEM following closely 
to the one in [11]. In the p- version of the FEM the polynomial degree p is used 
to control the accuracy of the solution while keeping the mesh fixed in contrast 
to the /;,- version or standard finite element method, where the polynomial degree 
is constant but the mesh size varies. Thus, the p-version is often referred to as 
the p-extension. The hp-method simply combines the h- and p- refinements. 

These different refinement strategies also imply different sets of unknowns or 
degrees of freedom: In the /i-version or the standard finite element method, the 
unknowns or degrees of freedom are associated with values at specified locations 
of the discretization of the computational domain, that is, the nodes of the 
mesh. In the p-method, the unknowns are coefficients of some polynomials that 
are associated with topological entities of the elements, nodes, sides, and the 
interior. 

For optimal /ip-convergence one should refine the mesh geometrically toward 
corners and let the degree of polynomial shape functions increase with distance 
from the corners. For an example of such a mesh see Figure |4j In the examples 
computed below, we have used a constant value of the order over the whole mesh. 

In the following one way to construct a p-type quadrilateral element is given. The 
construction of triangles follows similar lines. First of all, the choice of shape 
functions is not unique. We use the so-called hierarchic integrated Legendre 
shape functions. 

Legendre polynomials of degree n can be defined by a recursion formula 



{n + l)P„+i(x) - (2n + l)xPn{x) + nP„_i(a;) = 0, 
where Po{x) = 1 and Pi{x) = x. 

The derivatives can similarly be computed by using the recursion 
(1 - x^)P;^{x) = -nxPn{x) + nP„_i(x). 



The integrated Legendre polynomials are defined for x E [—1,1] as 




n = 2, 3, . . . , 



and can be rewritten as linear combinations of Legendre polynomials 



= {Pn{0 - Pn-2{0) , n = 2,3,.... 



The normalizing coefficients are chosen so that 
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Using these polynomials we can now define the shape functions for a quadrilateral 
reference element over the domain [—1,1] x [—1,1]. The shape functions are 
divided into three categories: nodal shape functions, side modes, and internal 
modes. 

There are four nodal shape functions. 



which taken alone define the standard four-node quadrilateral finite element. 
There are 4(p — 1) side modes associated with the sides of a quadrilateral (p > 2), 
with i — 2, . . . ,p, 



The internal shape functions are often referred to as bubble-functions. 

Note that some additional book-keeping is necessary. The Legendre polynomials 
have the property P„(— x) = (— l)"P„(x). This means that every edge must be 
globally parameterized in the same way in both elements to where it belongs. 
Otherwise unexpected cancellation in the degrees of freedom associated with the 
odd edge modes could occur. 

4.2. Solution of Lineeir Systems. Let us divide the degrees of freedom of 
a discretized quadrilateral into five sets, internal and boundary degrees of free- 
dom. The sets are denoted B.Dq, Di, iV°, and A^^, for internal, Dirichlet u = 0, 
Dirichlet u = 1, Neumann with Dirichlet m = in the conjugate problem, and 
Neumann with Dirichlet u — 1 in the conjugate problem, degrees of freedom, 
respectively. 

The discretized system is 



Niii, 77) = ^(1 - 0(1 - V), N,{i, r;) = ^(1 + 0(1 - V), 
N3i^,v) = ^(1+0(1 + ^), ^4(^,^) = ^(1-0(1 + ^), 




A 
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Taking the Dirichlet boundary conditions into account, we arrive at the following 
system of equations, using = 1, 

(Abb Asm ^bn" \ f / Abdi 

An^B A^ij^i Aj^fi^o j I X^vi I = ~ I ^AriDil 
A^OB AnOj^ji Ano^oJ \Xno/ \AnOb)^1^ 

For the conjugate problem, simply change the roles of Di ^ and Dq -h- N^. 
Note that Abb is, present in both systems and thus has to be factored only once. 

4.3. Evaluation of Contour Lines. Let u and v be solutions of the Dirichlet- 
Neumann problem and its conjugate problem, respectively. In computing the 
contour lines, the solutions and their gradients have to be evaluated at many 
points {x, y). Evaluation of the solution in /ip-FEM is more complicated than in 
the standard FEM. In a reference-element-based system such as ours, in order 
to evaluate the solution at point (a;, y) one must first find the enclosing element 
and then the local coordinates of the point on that element. Then every shape 
function has to be evaluated at the local coordinates of the point. This is outlined 



in detail in Algorithm 4.1 Similarly evaluation of the gradient of the solution 
requires two polynomial evaluations per one geometric search. 

Algorithm 4.1. (Evaluation ofu{x,y)) 

1. Find the enclosing element of{x,y). 

2. Find the local coordinates {^,f]) on the element. 

3. Evaluate the shape functions (f)i{C,,ri). 

4. Compute the linear combination of the shape functions Ci0i(^, 77), where 
Ci are the coefficients from the solution vector. 

Finding the images of the canonical domains is equivalent to finding the corre- 
sponding contour lines of u and v. Since both solutions have been computed on 
the same mesh, evaluating the solutions and their gradients at the same point is 



straightforward. In Algorithm 4.2 the two-level line search is described in detail. 



Algorithm 4.2. (Tracing of Contour Lines: u{x,y) = c = const.) 

1. Find the solutions u{x,y) and v{x,y). 

2. Set the step size a and the tolerance e. 

3. Choose the potential c. 

4. Search along the Neumann boundary for the point {x,y) such that u{x,y) = c. 

5. Take a step of length a along the contour line ofu{x,y) in the direction of 
Vv{x,y) to a new point {x,y). 

6. Correct the point {x, y) by searching in the orthogonal direction, i.e., Vm(x, y), 
until \u{x,y) — c\ < t is achieved. 



Conjugate Function Method for Numerical Conformal Mappings 



11 



7. Set (x, y) = (x, y) and repeat until the opposite Neumann boundary has 
been reached. 

4.4. On Computational Complexity. The solution time of a single problem 
can be divided into three parts, the setting up of the problem, the solution of 
the Dirichlet-Neumann and its conjugate problem, and the evaluation of the 
mappings. In short, in the absence of fully automatic mesh generators for this 
class of problems, the setting up of the problem remains the most time consuming 
part. We have implemented the algorithm using Mathematica 8. 
The time and memory requirements, in terms of degrees of freedom, have been 
reported for the Dirichlet-Neumann problems in In the examples below, 

the solution times vary from few seconds to at most two minutes on standard 
hardware (as defined in Table [T|. It should be noted that for comparable ac- 
curacy on polygonal domains, the Schwarz-Christoffel toolbox is superior to our 
implementation. 

Table 1. Effect of p and e on contour lines computations. In ge- 
ometry of Figure 10, ten contours of (radial) u{x,y) = Cu and (cir- 
cular) v{x,y) = c„ have been computed with a = 1/4. Times are 
normalized so that for p=8, e = 1/100, time = 1. The time units 
are thus Is and 46s for radial and circular contours, respectively. 
(Apple Mac Pro 2009 Edition 2.26 GHz, Mathematica 8.0.4) 



p 


e 


time for u{x, y) = Cu 


time for v{x, y) = Cy 


4 


1/100 


0.44 


0.43 


4 


1/1000 


0.41 


0.77 


4 


1/10000 


1.51 


1.19 


8 


1/100 


1 


1 


8 


1/1000 


1.00 


1.82 


8 


1/10000 


0.99 


2.66 


12 


1/100 


2.29 


2.31 


12 


1/1000 


2.26 


4.16 


12 


1/10000 


2.26 


6.07 



Estimating the computational complexity of the mappings is complicated, since 
in the end, the chosen resolution of the image is the dominant factor for the time 
required. In Table [T} the effects of the polynomial degree and the chosen tolerance 
on the overall execution time are described. As a test grid similar to one in 

Figure [10| has been computed by using a = 0.25 and e = 1/100, 1/1000, 1/10000, 
for p = 4,8, 12. Note that for the radial contours the effect of e is not as noticeable 
as for the circular ones due to contours and gradients being aligned. 
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5. Numerical Experiments 



Our numerical experiments are divided into three different categories. First we 
validate the algorithm against the results obtained by the Schwarz-Christoffel 
toolbox and the analytic formula. Then we study several examples of using our 
method to construct conformal mappings from simply (see Figures |6]-[9]) or doubly 



connected (see Figures 11-15) domains onto canonical domains, see Figure k5 



with the main results summarized in Tables [2] and |3| respectively. 



Table 2. Summary of the tests on simply connected domains. 



Accuracy is given as [log 



10 



M{Q)M{Q)\]. For the first two 



cases the moduli are known due to symmetry. 



Example 


ID 


M{Q) / M(Q) 


Accuracy 


Fi| 


2;ure 


Unit Disk 


5.1 


1/1 


-13 




6 




Flower 


5.2 


1/1 


-10 




7 




Circular quadrilateral 


5.3 


0.63058735108478 / 


-13 




8 








1.585823119159254 






Asteroid cusp 


5.4 


0.68435408764536 / 


-9 












1.4612318657235575 











Table 3. Summary of the tests on ring domains. Accuracy is 
given as [log^Q |1 — M((5)M((5)n, where the quadrilaterals are the 
cut domains. 



Example 


ID 


M{R) 


Accuracy 


Figure 


Disk in regular pentagon 


5.5 


See Table 


5 








10 




Cross in square 


5.6 


0.2862861647S 


8 


7473 


-9 




11 




Circle in square 


5.7 


0.9920378629010557 


-13 




12 




Flower in square 


5.8 


0.666955462334 


8065 


-8 




13 




Circle in L 


5.9 


1.0935085836560234 


-9 




14 




Droplet in square 


5.10 


0.897977509891 


8368 


-9 




15 





5.1. Setup of the Validation Test. Validation of the algorithm for the con- 
formal mapping will be carried out in two cases, first we compare our algorithm 
with SC Toolbox in a convex and a non-convex quadrilateral. In the second test 
we parameterized the modulus of a rectangle and map it onto the unit disk. 

The comparison to the SC Toolbox is carried out in the following quadrilat- 
erals: convex quadrilateral (fi; 0,1, 1.5 + 1.5i,i) and non-convex quadrilateral 
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Figure 3. Canonical domains and A^. on the left- and right- 
hand side, respectively. 

(VL; 0, 1, 0.3 + 0.3i, i), and line-segments joining the vertices as the boundary arcs. 
Then comparisons of the results obtained by the conjugate function method, 
presented in this paper, and SC Toolbox by Driscoll [8] are carried out. All SC 
Toolbox tests were carried with the settings precision = le-14. Comparison 
is done by using the following test function 

(3) ie^i{z) = \f{z)-g{z)l 

where / and g are obtained by the conjugate function method and SC Toolbox, 
respectively. The mesh setup of the quadrilaterals and the results of the test 
function ^ are shown in Figure |4] and |5| respectively. 




Figure 4. Geometric mesh of the convex (left-hand side) and the 
non-convex (right-hand side) quadrilateral used in computing the 
potential functions. 
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Figure 5. Comparison of the convex (left-hand side) and non- 
convex (right-hand side) quadrilateral between the conjugate func- 
tion method and SC Toolbox. Results are obtained by taking the 
logarithm (with base 10) of the test function (Isj). 



All our examples are carried out in the same fashion using the reciprocal identity 
Q and a quadrilateral Q. The test function is 

rec(Q) = |M(Q)M(g)-l|, 

which vanishes identically. See also [HI Section 4]. 

In the second validation test, we parameterized a rectangle with respect to the 
modulus M(Q) and map the rectangle onto the unit disk. The mapping is given 
by a composite mapping consisting of a Jacobi's elliptic sine function and a 
Mobius transformation. 

For every point {xj,yj) in the grid on the rectangle Rh, where Xj = j/lO and 
Vj — j = 0, 1, 2, ... , 10, we compute the error ||ej|| which is simply the 

Euclidean distance of the image of the initial point {xj, Uj) computed by the con- 
jugate function method and the analytical mapping. For a given modulus M((5) 
the values rec((5), max(||ej||), and mean(||ej||), where the latter two represent 
the maximal and the mean error over the grid are presented in Table |4] 

Note that our test function rec(Q) effectively measures the error in energy. Given 
the very high accuracy of the results obtained, we are confident that even though 
no a priori guarantees for pointwise convergence can be given, the second test is 
a valid indication of the global convergence behavior. 

5.2. Simply Connected Domains. In this section we consider a conformal 
mapping of a quadrilateral Q = {fl; Zi, Z2, z^, z^) with curved boundaries 71, 72, 73, 
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Table 4. The values of rec(Q), max(||ej||) and mean(||ej||) for a 
given M(Q). 



M(g) 


rec(Q) 


max( 






mean( 




) 


1 


8.08242 ■ 10- 


-14 


1.87409 


■ 10- 


-8 


5.56947 


10 


-10 


1.2 


6.35048 • 10- 


-14 


7.97889 


■ 10- 


-9 


7.49315 


10 


-10 


1.4 


5.52891 • 10- 


-14 


1.21851 


■ 10- 


-8 


6.90329 


10 


-10 


1.6 


8.85958 • 10- 


-14 


1.10001 


■ 10- 


-8 


7.90840 


10 


-10 


1.8 


9.72555 • 10- 


-14 


1.19005 


■ 10- 


-8 


7.31645 


10 


-10 


2 


9.41469 ■ 10" 


-14 


8.56068 


■ 10- 


-9 


7.67815 


10 


-10 



74 onto a rectangle Rh such that the vertices zi^ Z2, z^, z^ maps to 1 + ih, ih, 0, 1, 
respectively, and the boundary curves 71,72,73,74 maps onto the line segments 
7i; 72; 73; 74- We give some examples and applications with illustrations. Sim- 
ple examples of such domains are domains, where four or more points are con- 
nected with circular arcs. Some examples related to numerical methods and the 
Schwarz-Christoffel formula for such domains can be found in the literature, e.g.. 

Example 5.1 (Unit disk). Let Q be the unit disk. We consider a quadrilateral 
Q = {Q; zi, Z2, Z3, Zi), where Zj = e*^^ , 9j = (j — l)7r/2. Note that, because of the 
symmetry, it follows from ([T]) that the modulus is 1. The reciprocal error of the 
conformal mappings is 4.34 ■ 10"^*^. This example was first given by Schwarz in 
1869 [21]. 




Figure 6. Example of the conformal map of a square onto a disk, 
first obtained by Schwarz in 1869 |24j . 
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Example 5.2 (Flower). Let Q be the domain bounded by the curve 
(4) r(^) = 0.8 + tcos(n^), 

where < ^ < 27r, n = 6 and t = 0.1. We consider a quadrilateral Q = 
zi, Z2, Zs, Zi), where zj = r{9j), Oj = {j — l)vr/2; see Figure [7j As in Example 
|5.1 the modulus of Q is 1. The reciprocal error of the conformal mappings is 
3.74 ■ 10~^^. Several other examples of flower shaped quadrilaterals are given in 
[m Section 8.5]. 



Figure 7. Illustration of the flower domain and the visualization 
of the pre-image of the rectangular grid (Figure [3|. 

Example 5.3 (Circular Quadrilateral). In [TT] several experiments of circular 

quadrilaterals are given. Let us consider a quadrilateral whose sides are circular 

arcs of intersecting orthogonal circles, i.e., angles are tt/2. Let 0<a<b<c< 

271 and choose the points {1, e'", e*^, e^'^} on the unit circle. Let Qa stand for the 

domain which is obtained from the unit disk by cutting away regions bounded by 

the two orthogonal arcs with endpoints {1, e*"^} and {e*^, e*'^}, respectively. Then 

Qa determines a quadrilateral (Q^; e*", e***, e*'^, 1). Then for the triple (a,6, c) = 

(7r/12,177r/12,37r/2), the modulus M(Q^) = 0.630587351084775 and M{Qa) = 

1.5858231191592544. The reciprocal error of the conformal mapping is 1.68 • 
10-13 

Example 5.4 (Asteroid Cusp). Asteroid cusp is a domain Q given by a 

(5) Gc = {ix,y) : |x| < c, \y\ < c}, 

where c = 1 and the left-hand side vertical boundary line-segment is replaced by 
the following curve 

r{t) = -l + cos^t + isin^t, t e [-n/2,n/2]. 
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Figure 8. The quadrilateral {Qa; e*''/^^ ^'^rTr/u^ gi37r/2^ ^^^^ ^Yie 
visualization of the pre- image of the rectangular grid (Figure |3|. 

We consider a quadrilateral Q = {Q; 1 — i, 1 + i, + —1 — i). The reciprocal 
error of the conformal mappings is of the order 10~^°. The modulus M(Q) = 
0.68435408764536 and M(Q) = 1.4612318657235575. The domain is illustrated 
in Figure [9] 



Figure 9. Asteroid cusp domain with the pre- image of the rect- 
angular grid (Figure [3|. 

5.3. Ring Domains. In this section we shall give several examples of conformal 
mapping from a ring domain R onto an annulus A^. It is also possible to use the 
Schwarz-Christoffel method, see [13]. For symmetrical ring domains a conformal 
mapping can be obtained by using Schwarz' symmetries. 
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Example 5.5 (Disk in Regular Pentagon). Let Qhe a regular pentagon centered 
at the origin and having short radius (apothem) equal to 1 such that the corners 
of the pentagon are Zk = sec(7^/5)e2'^^''/^ k = 0,1,2,3,4. Let D(r) = {z E C : 
\z\ < r}. We consider a ring domain R = fi\D(r) and compute the modulus 
M{R) and the exponential of the modulus e'^^^^ Results are reported in Table 
[s] with the values e^*^^-* from Example 5] in the fourth column. 



Table 5. The values M(_R) and e 



M(_R) 



r 


M{R) 


exp(M(/2)) 


[3', Example 5] 


0.1 
0.4 
0.9 
0.99 
0.999 


2.35372035858745 
0.9674246001764809 
0.15070188000332954 
0.03276861064365647 
0.00934656029871744 


10.524652459913115 

2.631159438480101 

1.1626499971978235 

1.0333114143138304 

1.0093903757950962 


10.5246525 
2.631159439 
1.1626499972 
1.03331141431 
1.00939037579 




Figure 10. Disk in pentagon (r 
annular grid (Figure [3|. 



0.4) with the pre- image of the 



Example 5.6 (Cross in Square). Let Gab = {{x,y) : |a;| < a 
l^^l < b, \y\ < a}, and Gc as in (|5|, where a < c and b < c 
cross in square is a ring domain R = Gc\Gab 
error of the conformal mapping is of the order 10 
0.2862861647287473. 



see Figure 11 

10 



y\ <b}U {{x,y) : 
Then the domain 
The reciprocal 



The modulus M{R) 
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Example 5.7 (Circle in Square). Let Q be the unit disk. Then we consider a ring 
domain R = GA^, where c = 1.5, see Figure 12 The reciprocal error of the con- 

The modulus M(i?) 



Gc\^, where c 
formal mapping is of the order 10" 



0.9920378629010557. 




Figure 12. Disk in a square domain with the pre-image of the 
annular grid (Figure [3|. 

Example 5.8 (Flower in Square). Let Q he a. domain bounded by the curve 
Then we consider a ring domain R = Gc\^, where Gc is given by ^ and 



c = 1.5. See Figure 13 for the illustration. The reciprocal error of the conformal 



mapping is of the order 10"^. The modulus M{R) = 0.6669554623348065. 
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Figure 13. Flower in a square domain with the pre-image of the 
annular grid (Figure [3|. 

Example 5.9 (Circle in L). Let Li = {2; G C : < Re{z) < a, < lm{z) < b} 
and L2 = e C : < Re{z) < d, < lm{z) < c}, where < d < a, < b < c. 
Then L(a, b, c, d) = LiU L2 is called an L-domain. Suppose that B>{zq, r) = {z E 
C : 1^; — 1 < f'}- We consider a ring domain R = L{a,b,c,d)\3{zo,r), where 
(a, b, c, d) = (3, 1, 2, 1), zo = 8/5 + 2i/5, and r = 1/5. See Figure |T4| 

In order to better illustrate the details of the mapping, a non-uniform grid has 
been used. For the real component the points x are 

x = {k/10 : fc = 0,l,...,9}U {99/100,999/10000,9999/10000,1}. 

For the imaginary component the points y are chosen on purely aesthetic basis 
as: 

y = {k/10 : A; = 1,2,...,9} U 

{0.316225, 0.324008,0.327831, 0.329278, 0.331005, 0.687482}. 

The reciprocal error of the conformal mapping is of the order 10^^*^. The modulus 
M{R) = 1.0935085836560234. 

Example 5.10 (Droplet in Square). Let Q^i be bounded by a Bezier curve: 

r(t) = -i- U5f + 75t^ - 525t^ + 469) + (f - t e [-1, 1]. 

640 62 

Then the domain droplet in square is a ring domain R = Gc\Qd, where Gc in 

given in the first example concerning ring domains. For visualization, see Figure 



15 The reciprocal error of the conformal mapping is of the order 10 . The 



modulus M{R) = 0.8979775098918368. 
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Figure 15. Droplet in square with the pre-image of the annular 
grid (Figure [3]). 
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